Wave function for odd frequency superconductors. 
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We revisit the question of nature of odd-frequency superconductors, first proposed by Berezinskii 
in 1974. 1 We start with the notion that order parameter of odd- frequency superconductors can be 
thought of as a time derivative of the odd-time pairing operator. It leads to the notion of the 
composite boson condensate. 2 To elucidate the nature of broken symmetry state in odd-frequency 
superconductors, we consider a wave function that properly captures the coherent condensate of 
composite charge 2e bosons in an odd-frequency superconductor. We consider the Hamiltonian 
which describes the equal-time composite boson condensation as proposed earlier in Phys. Rev. B 52, 
1271 (1995). We propose a BCS-like wave function that describes a composite condensate comprised 
of a spin-0 Cooper pair and a spin-1 magnon excitation. We derive the quasiparticle dispersion, the 
self-consistent equation for the order parameter and the density of states. We show that the coherent 
wave function approach recovers all the known proposerties of odd-frequency superconductors: the 
quasi-particle excitations are gapless and the superconducting transition requires a critical coupling. 

PACS numbers: 



I. INTRODUCTION 

The discussion about possible symmetry types of 
superconducting order parameter A(k, r) (r denotes 
imaginary time) has drawn significant research interest. 
The conventional singlet (triplet) superconductor follows 
PTA(k,r) = A(k,r) ([PTA(k,r) = -A(k,r)]) un- 
der parity P and time T transformations. The singlet 
(PT = 1) and triplet (PT = —1) conditions can be satis- 
fied by either taking P = T = 1 and P = -1, T = 1 
for an even-in- frequency gap, or P = T = — 1 and 
P = 1, T = —1 for odd-frequency pairing. 

Although mainstream discussions of superconudctiv- 
ity are for even-frequency pairing, there is a growing in- 
terest in understanding odd-frequency pairing. The dis- 
cussion of unconventional pairing (P = 1,T = —1) was 
initiated by Bcrczcnskii 1 to explain the superfluid phase 
of 3 He. Although his proposal of triplet odd- frequency 
pairing could not explain the superfluid phase of 3 He, it 
certainly motivated a search of other possibilities of the 
pairing symmetries. Balatsky and Abrahams 3 later ex- 
tended the concept of odd- frequency pairing to the singlet 
superconductor (P = T = — 1). 

Although the realization of the odd-frequency pair- 
ing in current systems is still under debate, several re- 
ports consider this possibility in a a number of systems. 
Odd-frequency pairing in the Kondo lattice has been in- 
vestigated to study superconductivity in heavy-fermion 
compounds, 5 The proximity effects in a superconductor- 
ferromagnet structure 6 , a normal- metal/superconductor 
junction 7 and diffusive normal metal/unconventional su- 
perconductor interface 8 have been attributed to odd- 
frequency pairing. The p-wave singlet odd-frequency 
pairing is argued to be a viable pairing in the coexis- 
tence region of antifcrromagnctism and superconductiv- 



ity and/or near the quantum critical point in CeCu2Si2 
and CeRhIn.5. 9 In addition, hydrated Na^CoC^ is sug- 
gested to support an s-wave triplet odd- frequency gap 10 . 
Very recently, Kalas et al. 11 have argued that the boson- 
fermion cold atom mixture exhibits s-wave triplet odd- 
frequency pairing above some critical coupling at which 
the mixture phase separates. 

Motivated by the growing interest and possibilities 
of odd-frequency pairing, here we address the missing 
part of the odd-frequency superconductivity discussion: 
what is the wave function of the odd-frequency super- 
conductors? One might wonder how one can even ask 
this question given that superconducting correlations of 
two fermion operators in odd-frequency superconductor 
do not have an equal time expectation value? We as- 
sume (pretty safe assumption in fact) that any state, in- 
cluding odd-frequency superconductor, does has a many 
body wave function that captures superconducting cor- 
relations. Any state of matter has an associated wave 
function \tp) that captures the amplitude distribution of 
the particles forming this state. Hence we are asking ex- 
actly this question about the many body wave function of 
the odd frequency superconductors. Our wave function 
builds upon a long discussion 7 on the possible order pa- 
rameter and equal time composite operators that capture 
superconducting correlations of odd-frequency supercon- 
ductors in equal time domain. 

We propose a BCS-like pairing wave function for 
an odd-frequency superconductor, and study its conse- 
quences for the energy dispersion, superconducting order 
parameter, and density of states. The wave function, 
which describes a condensate of a spin-0 Cooper pair 
and a spin-1 magnon excitation, is consistent with the 
Hamiltonian suggested earlier in 2 to study odd-frequency 
superconductivity. We minimize this Hamiltonian with 
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respect to the proposed wave function and derive an ex- 
pression for the quasiparticle dispersion, a self consistent 
gap equation and the density of states. We find that a) 
the quasi-particle dispersion is gapless, b) the gap equa- 
tion has non zero solution only for a critical value of the 
coupling, c) the density of states is finite even for an en- 
ergy less than the gap energy, and d) the density of states 
is reduced at the gap edge compared to that of the BCS 
case. 

Before introducing the wave function and getting into 
the details of the minimization of the Hamiltonian, we 
would like to show that PT = 1 can be obtained by 
taking P = T = — 1 in 5 = singlet case. Any supercon- 
ducting order with translational invariancc, equilibrium 
and broken U(l) symmetry would result in an anomalous 
(Gor'kov) Green's function 

Fapfak) = (T T C ai k(T)c0,_k(O)}, (1) 

where a, (3 are spin indices. We assume that the transi- 
tion occurs only in a well defined representation. Thus, 
for S = singlet pairing, we may define 

F(r,k) = e Q(3 ^ Q/3 (T,k), (2) 

and for 5=1 triplet pairing, 

F(t, k) = (id^apTapir, k). (3) 

We now show the properties of F(t, k) under P and T 
transformations. For 5 = from Eq. (2), 

F(k,r) = e a p[0 r {c aik (T)cp _ k (O))-0_ T (c0 ,-k(0)c a ,kW>], 

(4) 

where 8 T is the Hcavisidc thcta function. 
We apply PT to this F: 

F(-k,-r) 

= e^[6 | -r{c M ,-k(-T)c l/i k(0)} - 6» r {c lyik (0)c M ,_ k (-T)}] 
= eju/[0-r(c M ,-k(O)c„,k(T)) - 6» T (c^k(r)c Al .- k (0))], 

(5) 

where in the last line we have used the fact that 
(TA(-t)B(0)) = (TA(O)B(t)) which agrees with the 
cyclicity of the trace, 

(A(-r)P(O)) = Ti{e~ HT Ae HT B) = 
Tr(Ae HT Be- Hr ) = (A(0)B(r)>. (6) 

Going back to Eq. 5, we permute fj, v, 
F(-K-t) 

= e |i „[e T <c > ,, k (T)c Wl -k(0)> - e_ T (c rk (0)c ftk (T))] (7) 
= P(k,r). 

All these properties of the Gor'kov function will be re- 
flected in the behavior of the gap function as well. There- 
fore, the gap function in general is even only under simul- 
taneous transformation: k — > — k (P) and r — > — r (T). 
We recall that PT = 1 is not only satisfied by P = +1, 
T = +1 but also by P = -1 and T = -1. The former de- 
scribes the BCS ,s-wave (even-frequency) pairing whereas 
the latter describes odd-frequency pairing. 



II. HAMILTONIAN AND WAVE FUNCTION 

When the idea of the odd-frequency pairing was first 
formulated for the singlet superconductor, an effective 
spin-independent interaction mediated by phonon was 
considered. 4 It was realized that this kind of interaction 
was unphysical for the singlet pairing. The problem was 
solved by considering spin dependent electron-electron 
interactions. Odd- frequency pairing posed another prob- 
lem related to the selection of the order parameter. In the 
BCS case the order parameter is generated from the ex- 
pectation value, F(r,t;r',t' — ► f) = (ip(r, t)ip(r', t)). But 
for the odd-frequency superconductor the equal-time gap 
vanishes since the gap is odd in frequency. This problem 
was solved by taking dF(r,t;r',t')/dt\ t _ t / as the equal- 
time order parameter. 2 

A Hamiltonian having a spin dependent electron- 
electron interaction was introduced by Abrahams et al. 2 . 
Using the equation of motion they derived an expression 
for dF(r,t;r' ,t')/dt\t-+t' ■ It was shown that the equal- 
time condensate for odd-frequency pairing is the expec- 
tation value of the product of a pair operator and a spin 
excitation operator. In what follows, we adopt this ap- 
proach, but for an odd frequency s-wave m = 1 triplet 
phase. We rewrite the Hamiltonian from Ref. 2 in the 
following form, 

H = £ kTC k T c kT + Y e kJ.4; c k| + Y ^S+S- 

k k q 

+ Y V kiqpC k +| T cl k +|^q" c -i+fi c i+fT' S 'p ' ( 8 ) 

klqp 

where ektl refers to the kinetic energy of the fl electrons 
measured from the Fermi energy, cj q is the magnon ki- 
netic energy, and Vki, qp is an attractive interaction which 
mediates the condensation, c kcr and c k( j creates and an- 
nihilates electrons at the state ker. 5 ± describe magnon 
excitations. Using this Hamiltonian, we propose a BCS- 
like wave function and study the superconducting state. 
The proposed wave function is written as 

= Yliukq + «kqC k+ | T ct k+ ^ j5+)|0), (9) 

kq 

where |0) represents the vacuum for both the electrons 
and the spin bosons. This wave function describes the 
superposition of the wave functions having two paired 
electrons with k+^ and — k+^ momentum and carrying 
opposite spins and condensed along with spin excitations 
(5+). i>k q («kq) represent the amplitude of the occupa- 
tion (or unoccupation) of these electron pairs with the 
spin excitation. 

There are key properties that explain this particular 
choice of variational function: i) 1-0} is a coherent state of 
composite bosons (c k+ a|C^ k+a ^5q )) that carry charge 

2e; ii) this wave function describes a coherent state that 
has broken U(l) symmetry associated with supercon- 
ducting condensate, as can be explicitly verified by using 
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c k — > cxp(i0)c k ; iii) Composite boson that condenses is 
not a simple Cooper pair 2 but contains two fermions and 
a spin-1 boson; iv) composite boson field has finite ex- 
pectation value in this state 



(V , |c k+ | T cl k+ q ; 5'+|V) = "kq^kq 



(10) 



and therefore \ip) is a mean field wave function for the 
composite condensate. 

The normalization of the wave function is given by, 

W#) = II(l U kq| 2 + kq| 2 <S-S + } q ) = 1, (11) 
kq 

which implies that |«kq| 2 + bkq| 2 (5~5 + ) q = 1 for all 
k,q. 

To make a next step we need to find the expectation 
value of the Hamiltonian (Eq. 8) with respect to the wave 
function (Eq. 9) and minimize it. Then we will proceed 
to derive the quasi-particle dispersion, density of states, 
and the self-consistent equation for the order parameter. 



We proceed by defining the two quantities A and E that 
will turn out to be the gap parameter and the energy of 
a composite excitation. 

Ak <i = ~\ E y kiq P sin20 lpV /<S-S+> qV /(S-S+> p , 

(16a) 



ip 



kq 



<S-S+)q)2+A kq 



(ek + i^ + ? ( ^ +)q)2+A ^ (16b) 



Then 



sin20 kq = 2u kq u kq W (S 5+} q 



A 



kq 



and 



=os20 kq = «5"5+) q - < = 



(17a) 



(17b) 



III. TOTAL ENERGY AND ITS MINIMIZATION where we have introduced the abbreviation 



The calculation of each term in Eq. 8 is shown in 
Appendix. Using Eqs. A2, B2, C2, Dl, the total energy 
can be written as 

E = 5> k+f +e k -H +c q < > <rs + ) q )|t- kq | 2 <s-s+) q 

kq 

+ V klqpWkqUkqflpttr p (5"5 + ) q (5"5 + )p. (12) 

klqp 

Following the BCS method, we choose Ukqj^kq such 
that they satisfy the normalization condition so that 
Ukq = sin# kq and i> kq = cos # kq / yj (5~5+) q . Then the 
expression for the energy reads 

E = J2 cos2 ^kq(e k+f + e k _ f + w q (5-5+) q ) + 

kq 

\ J2 y kiq P sin 20 kq sin 29 lp J (S-S+)qyJ (S~S+) P . 



aq= e k+| +e k _ f + ^ (g _^ (ig) 

Solving the normalization condition and Eq. 17, we can 
show that, 



1 



y kq 



2(5-5+), 



kq 



(19a) 



(19b) 



klqp 



(13) 



The minimization of the energy with respect to # kq gives 
dE 



d6 kc 



= - sin 20 kq (e k+ a + e k _a + w q (5 5 + } q ) 



^ V klqp cos 2(? kq sin 26 lp yj (S~ 5+ ) q yj (S~ 5+ > p = 0, 
ip 

(14) 

which can be rewritten as 



tan 26> kq 



Ei P *W sin 26 lpy /(S-S+} cl y/(S-S+) p 



BCS limit can be recovered at any stage of this analy- 
sis if we assume that spin correlators are factorized and 
have a peak at q = 0. This limit corresponds to the con- 
densation of spin field (5 _ 5 + } q = (5~) q (5 + ) q (5q,o. In 
this limit additional summation over q drops out and 
we recover standard BCS logarithm in sclfconsistcncy 
equation Eq.(16a), along with other features of BCS so- 
lution. This limit corresponds to the factorizitation of 
composite boson into product (tp\c^ +a ^c^_ k+ q,S^ \ip) — > 

^l4T ct -kilV'>W'l^l^}<5q,o. 



IV. ENERGY SPECTRUM 

Unlike the BCS case, i? kq is not a single-particle ex- 
citation energy. Therefore, we shall derive an expression 
for the energy required to excite an electron from the su- 
perconducting ground state. The excited state for an up 
spin is given by, 



e k+3 +£k-a +^ q (5 5 H 



(15) 



■0T 



I| (Wkq + W k qfc kq )]c k , + H T |0), 



(20) 



q : k#k' 
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where we have denned the composite creation operator 
^icq = c k+aT c -k+aj,'~'q ■ ^ e calculate the expectation 
value of the Hamiltonian Eq. (8) with respect to the ex- 
cited state wave function Eq. (20). The details are in 
Appendix E. The expectation value can be expressed as 

At 

(V-tI^t) = (^l#lV')+ek' + s + 1 ^-2Wk'q<S-S + ) q . 

" ^k'q 

(21) 

Using Eq. (19b). we can rewrite the above equation as, 



A£ T = e k ,+a 



Ck'q + ^k'q, 



(22) 



where AEf = (ijj^ \H\ipf) — (il)\H\ip) is the excitation en- 
ergy of the up spin electrons. AE^ can also be written 
as A£ T = k' • q/2m* - (w q /2)(5" S+) q + E k / q . Doing 
the same for the down spin excited state ipi, we find 
AE t = -k' • q/2m* - (w q /2)(S-S+) q + £ k ' q . 



V. DENSITY OF STATES 

The density of states (DOS) as a function of energy, 
N(E), is defined as, 

N ± (E) = J2 S[E (±^5 _ ^-{S-S+U + £ kq )], 



kq 



(23) 



where ± corresponds to up and down spins respectively. 
We numerically calculate the density of states for two 
cases of the magnon dispersion: 1) w q = q 2 /2M, and 2) 
w q = cu .. We set (S*-S*+) q = 1.0, and M = 10m*. The 
DOS for case 1) is shown in Fig. 1. 

In Fig la we show the DOS as a function of energy and 
order parameter. We have set a magnon momentum cut- 
off, q c = 0.25kp. We see that the DOS can be non-zero 
for energies less than the superconducting gap parameter; 
hence the DOS is gapless. The maximum of the density 
of state is always at the gap edge, but it is highly reduced 
at the gap edge compared to the BCS case. At E = 0, 
the DOS can be non-zero for small A. The calculation 
for a smaller q c (not shown in the figure) shows that the 
gap becomes more prominent in the DOS and spectral 
weight is transferred to the gap edge, similar to the BCS 
case. Hence q c — > reproduces the BCS results. In Fig. 
lb-d we have shown the plane cut of Fig. la for different 
values of A. For A = 0.1 (Fig. lb) we see that the DOS 
is non-zero for 0.05 < E < A. For A = 0.04 (Fig. lc) we 
see that the gap is completely closed and the excitations 
will be gapless. The effect is even bigger for A = 0.02. 

We also calculated the DOS using case 2): ojq = 
q 2 /2Af for q c > kp (the Fermi momentum) for a fixed 
value of A = 0.1. The result is shown in Fig. 2. In 
this figure we can see that the DOS almost closes the 
gap when q c = hp. As we increase q C7 the gap closes 
completely. Then the quasiparticle excitations become 



gapless. A still further increase in q c results in a finite 
DOS at E = 0. For q c > kp there is no enhancement of 
the spectral weight at the gap edge. 

The calculation of the DOS for w q = ujq also shows 
the similar density of state as discussed above for both 
q c = 0.25kp and q c > kp. 



VI. SUPERCONDUCTING GAP VS COUPLING 
CONSTANT 

The self-consistent gap equation (Eq. 16a) can be writ- 
ten as, 



A kq = 2 £ IT y/(SS+UyJ(S-S+) p , (24 ) 
i P lp 



where we have taken 



Vklpq = 



-V : |e k |<0.2^t 
: |e k | > 0.2^ 



(25) 



Then A kq = A = Ai p . The use of a more compli- 
cated interaction potential with a momentum depen- 
dence, would bring additional calculational complica- 
tions, which would not change the nature of the results. 



We denote p 



(wp/2)<5- 5+) p by f{p). We first 



perform the energy integral in Eq. 24 as follows, 



1 = 



i P 



(s-s+\ 



2 W(e P + /(p)) 2 



9 



= 9 



p 2 dp 



+■ A 2 
(S-S+) p de 



o v^p + f(p)) 2 + A 2 



J p 2 dp(S-S+) p log 



e c (p) + yj A 2 + e c {pf 
f(p) + x/A 2 + f{ P Y 



(26) 



where N(0) is the DOS in the normal state at the 
Fermi energy, g is the dimensionless coupling N(0)V/2ir 2 , 
e c(p) = ti^c + f(p) an d hjJ c = 0.2/x. If we assume spin 
correlator to have a sharp peak 8^ we recover BCS self- 
consistency equation from this equation. 

In the BCS c ase the gap equation is 1 = 
N(ti)V\og[{frw c + \JA 2 + hu; 2 )/A}. There is a solution 
for A for an arbitrary small value of A^(0)U due to log- 
arithmic divergence of the integral. In our case, in the 
presence of the magnon, the denominator will have some 
nonzero value because of the non-zero magnon energy. 
Then the right hand side can be made equal to 1 only for 
some critical value of g, as can be seen in the numerical 
evaluation discussed below. 



A. Case 1), lu p = p 2 /2M 

We solve Eq. 26 numerically for A as a function of the 
coupling strength g. We set m*/M = 0.1 and {S~S + ) P = 
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FIG. 1: The density of states (DOS) as a function of energy and superconducting order parameter in an odd-frequency 
superconductor. The DOS is normalized with respect to the DOS of the normal state. All the energies are normalized with 
Fermi energy of the system. The result in the upper panel is presented for the magnon momentum cutoff q c = 0.25. In this 
figure we show that the density of state is finite even for energy less than the gap energy. The maximum of the DOS is at the 
gap edge but the DOS is highly reduced compared to the BCS case. For smaller gap energy the DOS is completely gapless. 
Once the gap is closed the DOS starts to pile up at E = for smaller values of A. 
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FIG. 2: The DOS at fixed A = 0.1 as a function of energy 
for the magnon momentum cutoff q c > &f- The values of q c 
are 0.9, 1.1, 1.4, 1.7 times As we increase q c the gap in the 
density of states gradually closes up. For bigger q c the DOS 
piles up at E — 0. 



FIG. 3: The order parameter is numerically calculated for 
LOq — q 2 /2M. The magnon momentum cutoff is given by 
q c = Bk F , where B = 0.12,0.1,0.08,0.06, top to bottom. 
The order parameter is non zero only for critical value of 
the coupling g = N(0)V/2iv 2 . For a given value of N(0)V the 
larger A corresponds to the larger magnon momentum cutoff. 



B. Case 2), ui p — ujq 



1.0. The cutoff for the magnon momentum is given by 
q c = Bkp, where B varies between 0.12 to 0.06 in equal 
steps of 0.02. The result is shown in Fig. 3. In this figure 
we can see that a nonzero order parameter requires a 
critical coupling. 



The gap equation is again given by Eq. 26 but now 

fip) = p 2 /8m* - lo q (S-S + ) p /2. (27) 

We solve Eq. 26 numerically for A as a function of the 
coupling strength g. We fix the cutoff for the magnon 
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momentum to be O.lfci?. The result for various loq = Cfx 
where C = 0, 0.02, 0.04, 0.06, 0.08 is shown in Fig. 4. 
Again, the superconducting transition requires a critical 
coupling. 




0.4 



FIG. 4: The order parameter is calculated for to q = too- The 
magnon momentum cutoff is q c = 0.1fc_F. The result is shown 
for lo = Cfj, where C = 0, 0.02, 0.04, 0.06, 0.8, top to bottom. 
The order parameter is non-zero only for coupling exceeding 
a critical value. For a given value of N(0)V the bigger A 
corresponds to the smaller value of ojq. 



VII. MEISSNER EFFECT 



The Mcissner effect is one of the defining properties 
of a superconductor. The Meissner effect has been de- 
rived for the composite odd-frequency superconductor by 
Abrahams ct al. 2 . Here, we summarize the derivation 
given in that reference 

A superconductor shows the Meissner effect when the 
paramagnetic electrodynamic response is less than the 
diamagnetic response. The dc response is given by, 



ji(q) = -Qij(q)Aj(q), 

Ne 2 



Q?,-(q), 



(28) 



where A(q) is the Fourier transform of vector potential 
A(r), N is the electron density, and m is their mass. 
QiM) is given by, 



— e 
4m 2 



X]E k ^ / ^(^4( k +' r ) c 7( k -^)4(k-,o) C5 (k' + ,o)), 

kk' " 13 



(29) 



r 



where k± = k ± q/2. Q p can be evaluated near the crit- 
ical temperature T c by perturbation in the order param- 
eter A. The relevant Feynman diagrams of the current- 



J 



current correlation function for the Meissner effect are 
used. The analytical expression for q — > is 



<&(q)-QS(q) 



^[G 2 (k,Lu)G 2 (k\u;')-2G 3 (k,Lu)G(k\u;')}D(k + k\LU + uj'), (30) 
I 



.kk' 



where, G(k, u>) and D(k 1 lo) are the electron and magnon 
propagators. The condition for the Meissner effect is 
given by Q p — Q n > 0, which signifies the positive su- 
pcrfluid density in the superconductor. 

Situations with several models of the magnon propa- 
gators arc discussed. If the magnon propagator is mo- 
mentum independent, there is no contribution to (q) 
since the momentum summands are odd functions. So 
a momentum-dependent magnon propagator is used to 
discuss the Meissner effect. In the case of a static, spa- 
tially uniform magnon propagator having factorized form 
given by, D(q,i/) = —S^S^, the Meissner effect is found 
(Q p — Q n > 0). For spread-out (5-functions, the sign of 



Q p — Q n does not change, thus a positive supcrfiuid den- 
sity with a value between zero and the BCS value. Thus, 
it is shown that the composite odd-frequency supercon- 
ductors exhibit the Meissner effect. 



VIII. CONCLUSION 

In this paper we propose a BCS-like wave func- 
tion for the s-wave triplet odd frequency superconduc- 
tor. Our alternative approach to the odd-frequency su- 
perconductivity is based on the earlier discussion on 
composite bosons 2 . We present the wave function for 
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the odd frequency superconductor — Ilkq( w kq + 

t 'kqC ] J c+ a|C^ k+ q^S'+)|0}, , Eq.(9), that explicitly contains 

only the equal time operators and hence does not involve 
frequency or time domain. The wave function describes a 
condensate of a Cooper pair of spin S — and a magnon 
of spin S = 1. \tfj) does describe a coherent state that has 
nonzero expectation value for the composite boson op- 
erator, it captures the charge 2e condensate that breaks 
gauge symmetry and corresponds to the superconducting 
state. Naturally, since this \ip) describes odd-frequency 
superconductor, spatial parity P of this condensate is re- 
versed compared to the even frequency pairing opera- 
tors that corresponds to BCS condensate. Specifically, 
for the case we considered of spin triple S — 1 odd fre- 
quency condensate the spatial parity of the composite 
boson (c|(r)cj(r)5 + (r)) is P = +1 and hence this order 
parameter does posess all the quantum numbers inherent 
to the odd frequency S = 1 superconductor. 

We present a simplified model that captures the impor- 
tant features of the strong coupling theory developed for 
the odd-frequency superconductors and our results agree 
with the predictions of earlier studies: i) we show that 
the superconductivity requires a critical coupling. It was 
argued earlier 1,3 that a critical coupling is necessary in 



order to get the superconducting transition in the odd 
frequency superconductor, which we have also shown in 
this work, ii) we also derive the dispersion relation for the 
quasiparticles. We determine the density of states of the 
excitations. The density of states is very different from 
that of the BCS case. The gapless nature of quasiparticle 
excitations we find is also in agreement with earlier pre- 
dictions. The calculation of the density of states shows 
that it is always higher at the gap edge but its magnitude 
is highly reduced compared to the BCS case. For a range 
of parameters, unlike the BCS case, the DOS is finite for 
energies less than the gap energy and at E = it can be 
non-zero, hence odd-frequency supercoductor is gapless. 
We also argues how the BCS result is recovered by taking 
the magnon operator to condense and momentum cutoff 
q c = 0. 

Present discussion would be useful for the equal time 
formulation of the odd-frequency superconducting state 
and physical observables related to condensate. It also 
would be useful in elucidating the nature of condensate 
in odd-frequency supercondutors. 
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APPENDIX A: KINETIC ENERGY OF UP SPIN 
ELECTRONS 



It is convenient to rewrite X]k ek T c ikt Ck T as > 
(jJ Z) Z e kTC kT c kT = Z £ k+a T 4 c k+f T 



kq 



This is a trivial identity since we can shift k — » k + % 



(Al) 
a 

2 



and get the same result. 

We denote the (mn) component of the wave function 
as, \ip mn ) = (u mn + UmnC„ + a T cL m+ n_ L S'i l )|0). The ex- 
pectation value of the kinetic energy of the up spin elec- 
trons KE 1 = ^ Ekq(V'*|ek+a T 4 +aT c k+ | T |V') is given 
by, 



KE 1 = ^Z^k' g 'l £ k+a T cj c+ a T c k+ a T |i/'kq) 

kq 

= N X/k+f K'kq| 2 (S~S+} q . 



(A2) 



kq 



Here we use the normalization condition that 

Wn'^q'I^Wk.q) = <W <W ■ 

Then the kinetic energy of the up spin electrons is 

^T=Ekqek +f kq| 2 <S-S+)q. 
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APPENDIX B: KINETIC ENERGY OF DOWN 
SPIN ELECTRONS 

Using the same argument as discussed in appendix A, 
we rewrite, £ ki c kj, c kl as > 



q k 



kq 



(Bl) 

Then the expectation value of the kinetic 
energy of the down spin electrons KEy = 

W Ekq(^*l e -k+f iclk+s^k+f is given by 



KE i = ^I^kvl e -k+^ ct -k+s i C-k+ailV'kq> 



^E e -k + akq| 2 (5-5+) q . 



kq 



(B2) 



kq 



Then the kinetic energy of the down spin electrons is 
KE l = ^Ekq £ k-fkk q | 2 ^-5+) q . 



APPENDIX C: MAGNON ENERGY 

The expectation value of the magnon kinetic energy 
KE m = E q (V ; * l w q^q Sq IV') can be rewritten as 

KE m = (^)^(^ Y K5+S q |^ q ), (CI) 



k q 



which gives, 

KE m = Jr ^U^\V^\ 2 (S~ S + S~ S^q ( C2 ) 



kq 



APPENDIX D: INTERACTION ENERGY 

In the calculation of the expectation value of the in- 
teraction energy, Ei, it is easy to see that the product of 
only two states, kq and lp give non-zero contribution to 
the interaction term. All the other states are normalized 
to unity. Then, 



£/ = ^E<vhw>LApI^} 

klpq 

= E T/ klpq(V'kq^rpl4q 6 lpl^lpV'kq) (Dl) 

klqp 

r7 E v Uqp v kq u kqV\ P Ui p (S- S' + ) q (5" S+) p 



klqp 



where 6 tq = c fc+f T cL k+f ^5+ 



APPENDIX E: ENERGY OF EXCITED STATES 



The wave function of an excited state is, 

V~=[ II ("kq + WOR' + aTiO}- (El) 
q,k#k' 

Using the procedure of the Appendix A, we calculate 
the kinetic energy of the up spin electrons (KE^) with 
respect to the excited state wave function: 

^T = ^ E e k + akq| 2 (S-S + )q + e k ' + a (E2) 
q,k#k' 

where the restriction on k in the summation is inherited 
from the restriction imposed on the excited state wave 
function. £k'+a is due to the creation operator cL Ht 

2 K + 2 

which creates an up spin electron having unit probability 
of occupation in the state of momentum k' + If. We 
rewrite the Eq. E2 in the following form, 

KE \ = luE ek +f bkq| 2 ( 5 '~' S,+ )q + e k'+l - e k'+a|v k 'q| : 



kq 



(E3) 



Proceeding similarly, we show that the kinetic energy 
of the down spin electrons can be written as, 



KE ± = Jl E £ k-3 bkq| 2 (S S + )q - e k '-| K'c 
kq 



(E4) 



The kinetic energy of the magnon takes the following 
form, 



kq 



+ \2 

q' 



(E5) 



The interaction energy can be written as 

E I = E V klqpWl p ltr p U kqWfcq(5' _ 5' + )q('5 _ 5' + )p 

l (s-s + ) n (s-s + ) p . 



N 

klqp 

-2^V k 'l q] 
IP 



(E6) 



From Eq. 16a, we can show that 



A k 'q = - E Wipi\/(5^ + >,<S-5+) P . (E7) 
ip 

We use this relation in the right hand side of Eq. E6. 
The second term now gives +2u k , q u k 'q-\/ (S' _ S , +) q A k ' q , 
which, using Eq. 17a gives A k , q /£' k ' q . Then, 

1 A 2 , 

El = — E ^klqpWlpU 1 * p W kq Ukq(5' _ 5' + )q<5"5 + )p + 



klqp 



k'q 



(E8) 



Combining Eqs. (E3,E4,E5,E8), we get the result that 
we will use in DOS calculation, 



(ek'+3 + £ k'-§ +Uq(S S + )q)\v w ^\ 2 (S S + ) 



^k'q (E9) 
ql "5' + )q 



